arXiv:hep-th/051203 lv2 13 Dec 2005 


hep-th/0512031 


MCTP-05-100 


QCD recursion relations from the largest time equation 


Diana Vaman, York-Peng Yao 

Michigan Center for Theoretical Physics 
Randall Laboratory of Physics, The University of Michigan 
Ann Arbor, MI f 8109-1120 


Abstract 

We show how by reassembling the tree level gluon Feynman diagrams in a convenient 
gauge, space-cone, we can explicitly derive the BCFW recursion relations. Moreover, the 
proof of the gluon recursion relations hinges on an identity in momentum space which we 
show to be nothing but the Fourier transform of the largest time equation. Our approach 
lends itself to natural generalizations to include massive scalars and even fermions. 



1 Introduction 


QCD calculations are notoriously tedious if one is to follow the usual Feynman-Dyson ex¬ 
pansion in some commonly used, such as Feynman, gauge. Over the past few years, great 
strides have been made to simplify such endeavors. The results for the complete amplitudes at 
the tree or one- loop level can be quite compact. 

Following Witten’s proposal for a description of perturbative Yang-Mills gauge theory as 
a string theory on twistor space [ 1 ], and subsequent proposal for an alternative to the usual 
Feynman diagrams in terms of the so-called maximally helicity violating (MHV) vertices [2], 
a new set of methods was available for the computation of QCD amplitudes. The latest ad¬ 
vance in the form of recursion relations [3,4], in conjunction with the attendant rules for their 
construction, is particularly appealing. It is quite obvious from the flavor of such an approach 
that it bears on the cutting rules in held theory. In fact, some work at the one-loop level under 
the heading of cut-constructibility clearly points to the same origin [5]. These brief remarks 
certainly call for the possibility to develop the subject further in the context of quantum held 
theory and it is our intention to do so in this short article. 

As is well-known, unitarity of the S-matrix and the feasibility of an ordering of a sequence 
of space time points are intimately related. Indeed, the ordering need not be with respect to 
time, as is conventionally done. All that is essential in a perturbation series is that one must be 
able to separate the positive frequency and the negative frequency components in a propagator 
according to the signature of a certain linear combination of components Ax of the four vector 
between the two space-time points. For our purpose, a component rj ■ Ax of the light-cone 
variables will be a convenient start, where rj is a light-like vector. We shall rely on the existence 
of tubes of analyticity to continue such variables into the space cone, in order to incorporate a 
gauge condition for QCD. The resulting ordering is the equivalent of the largest time equations. 
We shall show that it is a consequence of these equations, when transcribed into momentum 
space, which give rise to a set of recursion relations. The outcome, for QCD in particular, 
is that one factorizes a physical amplitude into products of physical amplitudes, with some 
momenta shifted but still on-shell. This is the content of the BCFW recursion relations [3,4]: 

A(P, {Pi},Q , {Q,}) = £ A l (P, {PY) (p + y p) 2 A r (Q, {Qj}), 

where Ap , Ar are lower n-point functions obtained by isolating two reference gluons with shifted 
momenta, P = P — zr /, Q = Q + zr] with rj 2 = 77 • P = 77 ■ Q = 0, on the two sides of the cut. 
The shifting is necessary in order to preserve energy-momentum conservation. We would like 
to take this opportunity to point out that in so far as factorization is concerned, the masses 
of the internal propagators have no bearing. However, the demand that the shifted momenta, 
which will be called reference momenta, should be on-shell will force these external momenta 
to be light-like. 

We now turn to the important step of gauge fixing. In order to facilitate natural cancellation 
of terms at every level of a QCD calculation, the gauge that is most convenient for us is 
the space-cone gauge [ 6 ]. Here, the dependent degrees of freedom are completely eliminated 
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and only two hclicity components are left in the Lagrangian. To accomplish that, the now 
complexified null vector 77 is called upon. Loosely speaking, it is a direct product of two 
spinors. They are used on the one hand as the reference spinors separately for the two helicity 
components of the gauge field. On the other, they will be identified with two of the external 
momenta in a process. A further advantage of this gauge is that when we shift the momenta 
to obtain recursion relations, the dependence on momenta of the vertices will not be affected. 
Thereupon the factorization of the amplitudes is the same as that in a scalar theory. It is this 
special attribute which makes the program manageable. 

The plan of this article is as follows. In the next section, we write down the QCD Lagrangian 
in the space-cone gauge, where the auxiliary fields are eliminated. The number of diagrams 
contributing to a process will be drastically reduced. It is seen that the relevant propagator is 
basically scalar and that the vertices have good behavior under the kind of momentum shifts 
we shall make. This propagator will be decomposed into positive and negative frequencies in 
Section 3 according to light cone ordering. Sequencing space-time points in this ordering will 
be followed in Section 4 and the largest time equation will be summarized. 

To familiarize the reader with space-cone calculations and the underlying mechanism due to 
momentum shifting for factorization, we give several simple examples in Section 5. We want to 
emphasize that we obtain the recursion relations not only because the products of propagators 
satisfy certain algebraic identities, but just as importantly because the vertices also respect the 
kinematics in the shifted momenta. Furthermore, polarization factors have to be provided for 
the cut lines, so that the cut graphs are indeed physical amplitudes. The space-cone gauge 
fulfills all these demands. 

In Section 6 , we show that the propagator identity we used for five gluon amplitudes is a 
progeny of the largest time equation given earlier. It is this deep connection that we should be 
able to push the factorization program much further into the loop level. 

Section 7 is used to generalize the propagator identity to any number of space-time vertices 
which connect the flow of the two reference vectors. We shall show in fact that the propagators 
can have any masses, leading to a field theoretical proof of the recursion relations for massive 
charged scalars. At the same time this observation of course opens a new vista to include 
massive quarks, which we discuss in Section 8 . 


2 The space-cone gauge fixed Yang-Mills action 

Consider the four-dimensional Yang-Mills gauge theory, with the Lagrangian 

L = --Tr(d a A b - d b A a + i[A a , A b }) 2 . ( 2 . 1 ) 
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Following [ 6 ], we decompose the four-dimensional vector indices in a light-cone basis 

A a = (a, a, a + , a ~), ( 2 . 2 ) 

and we define the inner-product of two vectors by 

A ■ B = ab + ab — a + b~ — a~b + . (2.3) 
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Equivalently, we can choose to decompose a four-dimensional vector into bispinors: 

A * = = (“a a-) ■ ( 2A > 

The indices are raised and lowered using the northwest rule with the matrix 

Cap = Cap = -C°> = -C“i> = f° . (2.5) 

Null (light-like) vectors can be decomposed into a product of two commuting spinors 
(twistors): 

p«/3 = p y$ (2.6) 

Moreover, we can use the twistors to define a basis on the space of four-vectors 

P a = ( p\,Pa = \p),p a = [p\,Pa = \p \, (2.7) 


such that 

P = p + |+)[+| +p" |-)[-| +p|-)[+| +p|+)[-| • (2.8) 

With the normalization 

(H—) = [—f] = 1 (2.9) 

it follows that the components of a null four-vector P = \p)[p\ onto the twistor basis are given 
by 

p+ = (p-)[-p],p~ = (+p)[p+],p= (+p)[~p],p= (p—)[p+] - (2-10) 

As shown by [6], in a twistor formulation of the gauge theory a powerful simplification is 
achieved in the Feynman diagramatics by choosing the space-cone gauge: 


a = 0, 


( 2 . 11 ) 


followed by the elimination of the “auxiliary” component a from its equation of motion. The 
gauge fixed Lagrangian has now only two scalar degrees of freedom 


C = Tr 


d~ \ f d + 

-a 1 Da - — i ( — a + ) [a + , da~} —i[ — a~ 
2 V o J \ a 


[a da + ] + [a + , da } + [a + ,da ] —[a , <9a + ] 

a 2 


( 2 . 12 ) 

Choosing the space-cone gauge amounts to selecting two of the external momenta to be the 
reference null vectors for defining a twistor basis: |+)[+|, |—)[— |, such that the space-cone 
gauge fixing is equivalent to • A = 0, where the null vector N is equal to |+)[— |. The other 
ingredient which is needed in converting the essentially scalar Feynman diagrams arising from 
the gauge fixed Lagrangian (EH into definite hclicity gluon Feynman diagrams is inserting 
external line factors 

f+ = iW’ € “Fif (2 ' 13) 

for the positive, respectively negative helicity external gluons. The hclicities of the internal 
lines/virtual gluons are accounted for by the the scalar Lagrangian: a “ H —" helicity internal 
line corresponds to a a + a~ propagator, and vice versa. 
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3 The propagator 


For later purposes we explicitly construct a representation of the Feynman propagator, 
wherein a light-like four-vector is introduced as a parameter. Thus for a scalar held 1 one writes 


AO - y) 


1 r 1 iL .( X -y) 

i J (27t) 4 L 2 — ie 


(3.1) 


We will show that the propagator can be equally well represented as 

- (3.2) 

with y an arbitrary null vector. Appropriating the following common notations to the light-cone 
frame context: 


5+(p 2 ) = 5(p 2 )9(p + ), 5 0 2 ) = 5(p 2 )9(-p + ) (3.3) 

A + (x-y) = J -0^6 + ( pV'M, A-(x-!/)= f -0 ,F (?>"■ <*-»> , 

(3.4) 

we can rewrite the position space propagator El using a light-like ordering as 2 

A (x — y) = 9{(x — y) + )A + (x — y) + 9(—(x — y) + )A~(x — y ). (3.5) 

This can be seen, by first performing the contour integration over z 3 in (El 

/ z~ie e ^ ZV {X ~ y) = 2m6 ^ x ~ y 

r _dz_ e _ iz71<x _ y) = _ 2mQ ^ x _ , ( 3 . 6 ) 

J z \ 


under the assumption that the null vector y is equal to (y , 0,..., 0) and y >0, followed by 
an integration over p~. We end up with 


A(x - y) 


dp+ f dp i f 9(p+)9((x - y) + ) 9(-p + )9(-(x - y) + ) \ ip . (x _ v) 

iF/FoA-FI- + -FI- Y 


(3.7) 


1 Recall that in the gauge-fixed Yang-Mills action the gauge field was reduced to two propagating scalar 
degrees of freedom. So the scalar propagator constructed in this section is relevant also for gauge bosons. 

2 A more familiar form of the equation 1(1 511 . using a temporal ordering, is A(x — y) = 9((x — y)°)A + (a; — 
y) + 9{—(x — y)°) A~(x — y), where this time A ± (x — y) = J S(p 2 )9(±p°)e' lp ^ x ~ y \ 

3 Our conventions for the metric are mostly plus, and we define x ± = x Y* . 
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The latter expression can be reproduced starting from the usual Feynman propagator (EH), 
going to a light-cone frame and contour integrating over L . 

Note that the previous discussion can be extended to include massive particles. The minor 
extension involves replacing 5 ± (p 2 ) by S ± (p 2 + m 2 ), while rj remains a null vector. 

Finally, one can easily show that the propagator as defined in EH obeys the Klein-Gordon 
equation 

(— d 2 + m 2 )A(x — y) = -5 4 (x — y) . (3.8) 

i 

4 The causality (“largest time”) equations 

As stated in the Introduction, we shall show that the recursion relations are rooted in the 
largest time equation. To this end, we briefly revisit here the causality equations as derived by 
Veltman [7], but appropriately rewriting them in a light-cone frame. 

First, we introduce the following set of rules: 

-duplicate the Feynman diagram 2 N times, for N vertices, by adding circles around vertices 
in all possible ways; 

-each vertex can be circled or not; a circled vertex will bring a factor of i, and an uncircled 
vertex will bring a factor of (— i); 

-the propagator between two uncircled vertices is A (x — y), while the propagator between 
two circled vertices is the complex conjugate A*(x — y ); 

-the propagator between a circled Xk and an uncircled xi is A + (ay.—ay), while the propagator 
between an uncircled Xk and a circled xi is A~(ay, — xi). 

Clearly, the uncircled Feynman diagram is the usual one, while the fully circled diagram 
corresponds to its complex conjugate. 

The largest time equation states that the sum of all 2 N circled Feynman diagram vanishes: 

F{xi) + F*(xi) + F(ay) = 0 , (4.1) 

where F(ay) stands for the usual Feynman diagram, F*(ay) is its complex conjugate, and F(ay) 
is the sum of 2 N — 2 diagrams in which at least one vertex is circled and at least one is uncircled. 


Other causality equations can be obtained by singling out 2 vertices, Xk and ay. Let us 
assume x% < ay + . Then, one has 

6((xi - x k ) + )(F(xi) + F(k, x^) = 0 , (4.2) 

where F (k, ay) is the sum of all diagrams with k uncircled, but at least one other vertex circled. 
Similarly, one has 

0((x k - xi) + )(F( Xi ) + F(Z, Xi )) = 0 , (4.3) 

By adding these two equations one finds 

F(xi) = -F (k, l, Xi) - 6({xi - x k ) + )F(k, 1, ay) - 9((x k - ay) + )F(k, l, ay)), (4.4) 
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where F (k,l,Xi) is the sum of all diagrams with neither k,l circled, but at least one other 
vertex circled. F (k, 1, Xi ) is the sum of all amplitudes with k uncircled , but l circled and finally, 
F(k, l , Xi)) has k circled and l uncircled. 

By taking the Fourier transform of the real part of the position space Feynman diagrams 
- f dp(F(xi ) + F*(xi )) one obtains the imaginary part of the momentum space diagram. The 
cut graphs (Cutkosky rule) are derived by the use of ( 14 . 111 . The momentum space propagators 
between two vertices with both vertices uncircled is (14.111 — i/{jp 2 — ie), the complex conjugate 
expression if the two vertices are circled, and 2n5 + (p 2 ) if the momentum flows between an 
uncircled and a circled vertex. 


5 Reassembling Feynman diagrams into BCFW 
recursion relations 


The Feynman diagrams, as advertised, are the ones following from the space-cone gauge- 
fixed Lagrangian (ina . To set the stage for the recursion relations involving arbitrary tree 
level gluon n-point functions, we begin our investigation with the lowest ones. 

The 3-point function is the same as the 3-point vertex up to multiplication by the external 
lines polarization vectors. Otherwise, the 3-point function has only one peculiarity: to define 
it, one cannot select the 2 reference gluons out of the 3 external gluons. Take for instance 
(123) = (+ H—), with 1 selected as reference and an arbitrary null vector 4 = |4) [41. The 
3-point function is given by 


(+ H—) — e^~e 3 k 3 



(5.1) 


where the factor 1/(14) is inserted as a matter of normalization of the angular and square 
brackets. Using that 


(43) [23] — (41) [12], 


(5.2) 


the 3-point function acquires the standard googly-MHV expression 


(+ + -) 


[ 12] 3 

[23] [31] ’ 


(5.3) 


We now proceed to evaluate the 4-point functions. This is our first demonstration of factor¬ 
ization, which will be cast into a BCFW recursion. Let us consider (1234) = (4-h), with 

1 and 2 selected the reference gluons 


|1) = I”), |1] = H (5-4) 

|2) = |+), 12] = |+] . (5.5) 
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In this basis, the polarization vectors are 


+ = \-pL = M 

P (+p) (2 p) 

_ _{+p) _ (2 p) 

6p [-P] [1 P] ’ 

and the 4 components of the momenta decomposed in the bispinor basis read 

P + = [~P\(P-) = [ l P\{P l ) 

P~ = [+p\(p+) = IM(p 2 ) 

P = ~[-p](p+) = ~[lp](p2) 

P = -[+P\{P~) = *[2p]](pl) • 

Notice that the only non-vanishing components of the reference glnon momenta are 

pt = i, pi = i • 


(5.6) 

(5.7) 


(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 


We now pose the question what happens with the Feynman diagrams if for later purposes we 
decompose the //-shifted reference gluon momenta, with the choice 

7 — |2> [1| , (5.13) 


in the same twisor basis ? Concretely, we define the //-shifted reference gluon momenta to be 


Pi = Pi + zr), P 2 = P 2 


zrp 


(5.14) 


or, in components, 

Pi = 1, Pi = -*[21]<21> (5.15) 

Pt = 1, P 2 = z[21}(21). (5.16) 

Thus the only changes to the reference momenta enter through the component p. This is 
particularly important, since the vertices in the space-cone gauge turn out to be independent 
of p, as it can be seen by inspecting the Lagrangian (HUD - 

We conclude that any shift of the external momenta as in fl 5.14 ) is of no consequence for 
the vertex-dependence of any Feynman diagram, and leaves an imprint only over the internal 
line propagators. 

The 4-point function in the space-cone gauge is given by a single Feynman diagram 


which is represented by the left hand side of Fig.l. With the choices already made in terms of 
reference momenta 4 


(1234) = (+ - -+) 


- + 2 
e 3 e 4 PK 


1 1 

\P}> ( 12 ] ( 12 ) 


= iKl IN [131 (32) [141 (42)_i_i_ 

(13] (24) 11 { ' 11 ( ' [14] (41) [12] (12) 
(23) 2 [14| 

_ (41) (12) [12] ' 

4 The factor of 1/2 in the propagator is due to a peculiar normalization of [6]. 


(5.17) 

(5.18) 

(5.19) 
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Figure 1: Factorization of the 4-point function 

Using that [12] (23) = —[14] (43), which follows from momentum conservation one recovers the 
simple MHV expression of the 4-point function 

d 234 ) = <+ ~ -+) = (12) (34)( 41) ■ < 5 - 20 > 

On the other hand, we may choose to split the 4-point function into lower on-shell ampli¬ 
tudes, 3-point functions, as indicated on the right hand side of Fig.l. The internal leg K has 
been put on-shell by a shift of the two external reference gluons as in (EH1, where 

- = P » ( I4 >M < 14 > (52 n 

2 rrP u (24) [14] (24)' 1 

Notice that because there is no difference between the 3-point function and 3-point vertex, 
other than the multiplication by external line polarizations, and because we have essentially a 
scalar field theory, one can insert freely factors (=1)- Recalling that as argued before, the 
shift (15.141) in the external momenta does not modify the vertices, we see that the factorization 
into 3-point functions according to Fig.l. is trivially realized. Furthermore, if we follow the 
same steps we made earlier in converting the 3-point functions into MHV vertices, we arrive at 
the BCFW result. 

The 5-point, function is the first non-trivial example in which we invoke an identity rooted 
in the largest time equation. To begin with, it is given by the sum of three Feynman diagrams 
which are represented on the left hand side of Fig.2. The right hand side contains terms corre¬ 
sponding to having cut the 5-point function in all possible ways such that the reference gluons 
1,5 are on opposite sides of the cut. Moreover, we shift the reference external gluons such that 
the cut line is on-shell. Thus the cut diagrams become on-shell amplitudes. Moreover, the cut 
diagrams are multiplied by the propagator of the line that was cut. 


We will establish an identity relating the 5-point function, as computed from Feynman dia¬ 
grams, to the tree amplitudes of 3 and 4-point functions as indicated by Fig.2. This, of course, 
is nothing but the statement of the BCFW recursion relation, applied to this particular 5-point 
function. 

In particular we will show that the top Feynman diagram is the sum of two such cut diagrams 
(A+B), and that the middle and bottom Feynman diagrams are equal, respectively, with C 
and D. The equality of the two sides of the last two diagrams is obvious from the fact that, as 
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Figure 2: Factorization of the 5-point function 


argued before, the vertices on the left and right hand side of Fig.2 are the same, irrespective of 
having shifted the external reference momenta in order to put the cut line on shell. Moreover, 
the propagators on the right and left hand side of the bottom two diagrams coincide as well. 

The attentive reader could observe that B+D add up to zero. After giving an obvious 
common factor, B+D add up to the (+ + -I—) 4-point function. The reason why we have to 
consider two Feynman diagrams to recover the 4-point function, as opposed to our previous 
calculation, is that we have selected only one of the four external gluon momenta as reference 
vector. This means that we have to consider both the s and t channel. Nonetheless, the sum of 
these two channels is zero, as it corresponds to having all but one external gluons of the same 
helicity. Thus indeed, the BCFW recursion relation amounts to including only the terms A 

and C, corresponding to a factorization of the (+ + -f-) amplitude into (+ H—)(+ H-) 

amplitudes. 

In the derivation of the recursion relation directly from Feynman diagrams, it is useful to 
keep all possible terms that arise from cutting an internal line of all Feynman diagrams, such 
that the reference external gluons are on opposite sides of the cut. To show that the top 
Feynman diagram equals A+B, once we factored out the vertices (which are the same on the 
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left and right side of Fig. 2), amounts to proving the following identity between propagators: 


1 1 
p2 p2 

-*12 -*45 


p n n s 


1 1 

pfpl' 
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where we defined the shifted reference gluon external momenta 


P\ = Pi + zr,), 
P{ = Pi + zr), 


P§ = P 5 - zr) 
P\ = Pb-zr). 


z, z are such that we put the internal lines K , L , respectively, on-shell 

( 12 ) 


* = 


z = 


p 2 
-'12 


2 v ■ P 12 (52) 

(45) 


p 2 
- r 45 


(5.22) 


2r? • P 45 (15) ’ 

and where the null vector 77 is defined, as before, with respect to the reference gluons 

r7 = |+) [—| = 15) [1| . 

Despite the fact that we can prove 1ET221) by going into the bispinor basis, we find it 
simpler and easily allowing for generalizations to stay in momentum space. Using that 


p l 2 = -2 zr) ■ P 12 , P 45 = 2+77 • P 


45 


P h = ' Pl2 ’ P 4 S = 2 ( £ ~ *)r) ■ P 45 , 


(5.23) 

(5.24) 

(5.25) 

(5.26) 
much 

(5.27) 

(5.28) 


il5~22l) becomes a trivial algebraic identity 


1 

zz 




(5.29) 


This completes the proof of the BCFW recursion relation from Feynman diagrams for the 
5-point function and highlights the pattern that we will encounter for an arbitrary n-point 
function. 


6 The recursion relations and the largest time equation 

There is yet another way to address the identity (15.2211 by rewriting the propagators in 
position space and next recognizing the Fourier transform of the (i.e. momentum space) largest 
time equation (S3). Reinstating the usual ie prescription in the momentum space propagators, 
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and multiplying 1ET221I with the total momentum conservation 5-function, the right-hand-side 
of (15.221) becomes: 


S(Pi + --- + P 5 ) 1 

p2 p2 

- r 12 -*45 


dL i(Pi+P2+K) Xl i(P 3 -K+L)x 2 UPj+Pn-pxa 
2vr) 4 

1 1 

K 2 — ie L 2 — ie 

e I - *»>*<*» - »(a,) 


X\dx 2 dx 3 


■ dK 


The shifted propagators which appear on the left-hand-side of can be cast into 

5(Pi + --- + P 5 ) 1 


p2 

- r 12 


P 2 

45 


dx\dx 2 dx 3 


[ dz f dL c i(P 1 +P 2 +k)x lc i(P 3 -k+L+zr l )x2 c i(P4+P5-L-zv)x 3 


(2vr ) 4 J J (2vr) 4 

1 f 5~ {{k — zrj) 2 ) 5 + {{k — zr ]) 2 ) 


= i 


L 2 — ie\ z — ie z + ie 

2 1 dx 1 dx 2 dx 3 e i{Pl+P2+)xi+iP3X2+i ^ +P5)x3 


(e((x 1 -x 3 ) + )A + (x 1 -X 2 ) + 0((x 3 -x 1 ) + )A (xi-x 2 )Ja(x 2 -x 3 ), (6.2) 


using the same z-parametrization which we have introduced in Section 3. For concreteness we 
have chosen, as before, rj = ( r )~, 0,0) with r]~ > 0. Clearly, by integrating out first xi,x 2 , x 3 and 
using the delta-function to integrate over k, L , followed by a ^-integration using the remaining 
delta-function 5((Pi + P 2 — zr]) 2 ), we recover the left-hand-side of (16.21) . On the other hand, if 
we choose to shift the integration variable from k to k — zr /, and we next integrate over z, then 
we find the result given in the last line of (E3- 

Similarly, the other term on the left-hand-side of can be written as 


S( Pl + ... + p 5 ) 1 


P 2 

12 


p 2 

- r 45 


dx\dx2dx 3 


dK f dl 


(2tt) 4 ,/ (2tt) 4 
1 / 5~ ((/ — zrj) 2 ) 5 + {(l — zrj) 2 ) 


e 'i(Pl+P2+K+zr])x 1 e i(P 3 -K-zri+l)x2 e i(P4+P 5 -l)x 3 


K 2 — ie \ z — ie z + ie 

= i 2 I dxidx 2 dx 3 e i( ' Pl+p ^ xl+iP3X2+i< ' Pi+p ^ X3 


^6»((xi - x 3 ) + )A + (x 2 - x 3 ) + 0((x 3 - xi) + )A (x 2 - x 3 )^A(x 1 - x 2 ). (6.3) 
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Thus the identity dE 22 t becomes 


dxidx2dxse i< ' Pl+P2 ' >xl+iP3X2+ii ' P4 ‘ +P5 ' )X3 


A(aq - x 2 )A(x 2 - x 3 ) 


- (^9((x 1 -x 3 ) + )A + (x 1 -x 2 ) + 9((x 3 -x 1 ) + )A (x 1 - x 2 )^j A(x 2 - x 3 ) 
Q((x i - ^ 3 ) + ) A + (x 2 - x 3 ) + 0((x 3 - x\) + )A~(x 2 - x 3 )) A(xi - x 2 ) 


= 0. (6.4) 


There is one more step that is needed in order to show the relationship between (FT221) and 
the largest time equation. From (@3D, with xi, x 3 the two vertices that are singled out, we have 

A(x ± - x 2 )A(x 2 - x 3 ) = A~(xi - x 2 )A + (x 2 - x 3 ) 

+ 9((x i - x 3 ) + ) ^A + (xi - x 2 )A(x 2 - x 3 ) - A*(xi - a: 2 )A + (a:2 - x 3 )^ 


+ d({x 3 -Xi) + )( A(xi -x 2 )A (x 2 — x 3 ) — A (x x - x 2 )A*(x 2 -x 3 )j. 

(6.5) 


To show how (B is related to (16.51) we first rearrange the right-hand-side of (16.51) using 

A*(x -y) = 9((x - y) + )A~(x - y) + 9{(y - x) + )A + (x - y) , ( 6 . 6 ) 

such that it becomes equal to the Fourier transform of the right-hand-side of equation (D- up 
to the following two extra terms: 9((xi~x 3 ) + ) A + (xi—x 2 )A + (x 2 —x 3 ) and 9{{x 3 — Xi) + )A _ (a;i — 
x 2 )A~(x 2 — x 3 ). These terms in fact are zero as the product of the three distributions has zero 
support. The easiest to see this is to evaluate the Fourier transform 

J dxidx 2 dx 3 e l ( Pl+P2 ' >Xl+ ' lP3X2+l( ' P4+P5 ' >X:i 9((xi — x 3 ) + )A + (a:i — x 2 )A + (x 2 — x 3 ) (6.7) 

by rewriting the step function as a ^-integral, followed by the integration over x±,x 2 ,x 3 , to 
arrive at 


S(Pi + ...P 5 ) f dz—^—S+ttP! + P 2 + zrj) 2 )5+((P 4 + P 5 - zyf) . (6.8) 

J z-ie 

It is clear that no z can satisfy the simultaneously the two delta-function constraints. This 
completes the proof that the algebraic identity which was found by reassembling the Feynman 
diagrams into the BCFW recursion relations arises from the more fundamental largest time 
equation. 

Before closing this section, we need to add some remarks to explain some more what we 
have accomplished. The discussion above of the largest time equation is taken for real p, since 
it is only for such values that we know how to order a sequence of space-time points in r/ • x % . 
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Now that we have the largest time equation, let us Fourier- transform it into momentum 
space, appropriate for a physical process under consideration and keeping 77 as a variable. Then 
we obtain a set of shifted momenta, as in (15.231) . We then analytically complexify 77 . For tree 
level, this is certainly possible and justifiable, because the dependence on it is only in some 
algebraic functions of propagators. At the loop level, we need to invoke the analysis of axiomatic 
field theorists [ 8 ], which states that there are tubes of analyticity to allow this extension and 
to lead to complexified unitarity relations. We now identify these propagators with the ones 
which we need in the space-cone gauge to carry on with the analysis. 


7 The general case 


To exploit the full generality of the problem, we derive an identity satisfied by the momentum 
space scalar propagators working under the assumption that we deal with massive propagators, 
with arbitrary masses. 

Consider a graph at the tree level with n! vertices and m! external lines which are on-shell. 
As our convention, we take them to be all outgoing. We single out two of these lines which do 
not land on the same vertex as reference vectors and call them p a and pb . For a tree graph, 
there is a unique path through some of the internal lines which connects p a to pb . We shall 
denote the vertex at which p a emanates aq, and that for pb x n in their space-time labels. The 
vertices in between are aq, % — 2,... n — 1. There are then n vertices in this segment of the graph 
and therefore n-1 internal lines. Our consideration for the time being will be on this segment. 
The internal lines carry momenta < 77 , joining aq to aq+i ,77 — 1 > i > 1. How and what other 
lines enter or leave these vertices need not concern us at this point. For each 97 , we associate a 
propagator 


A(aq - x i+ i,mi) 


d A qi 


0 iqi-(xi-x i+ 1 ) 


ie 


J (2n)Hqf + mj 
0(77 • (aq - aq+i))A + (aq - Xi + i,rrii) + 6>(—77 • (aq - x i+1 ))A ~(aq 


1 1 rrii J 


(7.1) 


where 


A + (aq - aq+i, rrii) 
A “(aq - Xj+i, rrii) 


d A qj 

(2n)H 

d 4 qi 

(27t) 4 7 


e^< Xi - Xi+l) 8 + (^+ m 2 i), 
e m-( x i~ x i+i)§- _|_ rnf) , 


(7.2) 


and where 77 is a light-like vector. Due to momentum conservation, each 37 is expressible in 
terms of external momenta, and in particular it has a component +p ai or equivalently —pb- 
The factorization procedure is to cut these g* successively by shifting them by zr/. The 
on-shell conditions 


Qi + = 0 , qi = qi- z 


(7.3) 
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will give us a set of solutions, points in the complex plane, namely 


Q? + rnj 

2 ? 7 ' Qi 


(7.4) 


for each n — 1 > i > 1. More precisely stated, the factorization amounts to splicing the graph 
into a sum of products of two on-shell graphs with shifted momenta {p a — Ziq,...,qi} and 
{—qi, ■ ■ ■ ,pb + ^ 77 }, where ■ • • stand for the other momenta in the left graph segment and 
similarly for those in the right graph segment, with the propagator as the partition. The 

reason that p a and pb are shifted is because we need to conserve the overall momenta on the 
left and the right segment separately to make them into physical amplitudes. We must demand 
on-shell conditions for the shifted p a q, with the same masses, which give 


(Pa ~ ZiV) 2 + ml = 0, 


(Pb + Zip) 2 + m 2 b = 0, 


or 

Pa ■ rj = 0, p b ■ q = 0. 

As q is light-like, these conditions clearly do not allow p aj b to be time-like. Therefore, the two 
reference vectors must also be light-like, m a = mb = 0 . 

The identity which we want to establish is 


1 1 
q\ + m\ql + 


1 _ 1 1 1 
gLi+ m n-i qi + m\ (q 2 - Zxq ) 2 + ml (g„_i - z t q ) 2 + m 2 ^ 

11 1 
{qi ~ z 2 q ) 2 + m{q% + m\ (<?„_ 1 - z 2 q ) 2 + m 2 n _ x 

+ . 

1 11 

(Qi - z n -iq ) 2 + m\ (q n - 2 - z n _ 177) 2 + m 2 n _ 2 q 2 n _ x + m^_ 1 ' 

(7.5) 


The proof of this identity is quite simple. For n — 1 >i ^ j > 1, we write 

(ft - zjp) 2 + m 2 = q 2 + m 2 - 2 Zj q ■ qi. (7.6) 

Then, using the on-shell conditions for the shifted internal momenta, which is tantamount to 
making cuts, we have 

Qi + m 2 = 0 -> q 2 + m 2 = 2z t r] ■ q t . (7.7) 

Together, they yield 

(ft - Zjq) 2 + m 2 = 2q ■ qi(zi - Zj) . (7.8) 
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Putting these together, we see the identity holds if one can show 


(~l) n 

Z1Z2 ' • • Zn-j. 


1 

Zl(Zi - Z 2 )(z 1 -z 3 )--- (Zi - Z n _ i) 

1 

(Z2 - Z 1 )z 2 (z 2 - z 3 ) • • • (z 2 - z n - 1 ) 


(Zn —1 Z\)(z n —i Z 2 ) • • • 2)-2'n—1 


(7.9) 


This is so, because dZ3 is just a formula of partial fractioning, or it is just a statement that 
the integral 

f _ ± _= 0 

J z(z - Zi){Z - Z 2 ) ■ ■ ■ (z — Zn-!) 

for a complex variable z over a contour which encloses all the poles. 

Notice that eqn. m is precisely the identity needed to reassemble a generic tree level gluon 
Feynman diagram into lower on-shell amplitudes, as shown in Section 5. The reason for this is 
that, as we argued before, the vertices which enter the Feynman diagram and the corresponding 
lower n-point functions are the same , being insensitive to the shift of the reference gluons. Also 
the external line factors that have to be inserted on the cut lines to recover the lower on-shell 
amplitudes cancel pairwise, i.e. e + e _ = 1. Then, one is left to prove only an identity involving 
the momentum space scalar propagators. This is the same as Q. with all propagators being 
massless rrij = 0. The combinatorics work out properly to reproduce the BCFW recursions. 
We have also checked these points explicitly for the six gluon amplitudes. 

Furthermore, the arguments presented in Section 6, relating the momentum space identity 
m to the Fourier transform of the corresponding largest time equation @3D, with the sin¬ 
gled out vertices corresponding to those of the external reference gluons, can be easily carried 
through. 

This completes our purely field theoretical proof of the BCFW recursion relations. In the 
process, we have identified the underlying principle behind them in the form of the largest time 
equation. 


8 Adding massive scalars and fermions 

Establishing recursion relations to include charged massive scalars is straightforward in 
our framework. The current interaction term ( < F*(9 /J f h — Qd fl Q*)A fl has only d + ,d~ and d 
derivatives, without the dangerous 3 which would have been sensitive to the shift of the external 
momenta, because of the space-cone gauge (a = 0). The quartic interaction A^A 11 has no 
momentum dependence. Therefore the vertices are unchanged under the zg shifts. Besides 
there is no external line factor for scalars. This implies that we have met all the requirements 
to accommodate recursion relations as stressed over and again. 


15 








A more interesting and natural generalization is the inclusion of fermions in the tree level 
recursion relations. Consider the Lagrangian of minimally coupled massive fermions 

£f = ^ - mj)^i , (8.1) 

i 


where V M is the gauge covariant derivative and the Dirac matrices obey the usual anticom¬ 
mutation relations {7^, 7^} = — 2 g fiu . Then the recursion relations which are formulated with 
the two reference gluons connected by a path which includes a fermionic line are based on the 
following identity involving momentum space fermion propagators 5 


_... y ^-2 _ 

ii + m 1 j 2 + rn 2 j n -i + m n _i 

= m i~^i + 1 7 m 2 ... 7 ^-2 _ 

(qf + mf) fa ~ zih + m 2 j^i-i 

+ 1 + _ 

ii - z 2 fi + mi (ql + ml) ^ n -i 


1 

£ 1 ?) + m n -i 

1 

£2 fl + m n - 1 


^1 - Zn-ifj + mi 


-y/U . . . -yMn-3. 


in-2 


^n-iW + m n - 2 


-7 


l^n — 2 


mn -1 - jn-1 + Z n -.\ fj 


( 8 . 2 ) 


A Dirac matrix inserted between two fermion propagators corresponds to a cubic interaction 
vertex with a gluon held. For the case when this gluon held is an external one, one must insert 
external line factors and thus contract the space-time index of the Dirac matrix with that of 
the corresponding polarization vector. If the gluon held corresponds to an internal line, then 
we must contract the space-time index of the Dirac matrix with that of the gluon propagator. 
At this stage it is important to stress that we are in the space-cone gauge, such that rj ■ A = 0, 
which means that jfl. = —/l//. Thus in ( 17 . 51 ) we consider only insertions of Dirac matrices which 
effectively anticommute with the null vector 77. 

The proof of ( 18 . 21 ) relies partly on one identity which we have already established, namely 
m- First rewrite the fermionic propagators such that all denominators will correspond to 
scalar propagators. Next cancel out in the numerator all terms with at least one insertion of 
ij. This can be done, since ijfi = 0 , corroborated with the previous observation that 1/] and 7^ 
anticommute. Finally, we need to employ one other algebraic identity, namely 


dz 


(z - zi)(z -z 2 )...(z- Z n —\) 


= 0, for m < n — 2 


(8.3) 


where the integral is evaluated over a contour which encircles all the poles. As mentioned many 
times before, the structure of the vertices is unchanged by the shift with zjj in the momenta 
of the reference gluons. We complete the proof of the recursion relations by observing that 
each term in (Q corresponds to a factorization into lower on-shell amplitudes, such that each 

5 For concreteness we considered here a quark line. Otherwise, for an antiquark the signs in fj — Zjj] + rrii 
must be flipped to —ji + Zjf) + m*. 
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propagator belonging on the path that connects the two reference external gluons is cut and 
put on-shell, accompanied by the corresponding shift of the external gluons. The left and right 
amplitudes are multiplied by the propagator of the line which was cut. The reason why in the 
recursion relations the “cut” fermionic lines have an extra factor rn t — j t + Zij] has to do with 
the fact that this corresponds exactly to the appropriate insertion of the external line factors. 

Thus we have provided a field theoretical proof of the recursion relations for massive scalars 
and fermions found by [9,10]. 


9 Conclusions 

To summarize our results, we have shown that the BCFW recursion relations can be proven 
starting from the standard gauge theory Feynman diagrams. Perhaps the only less familiar 
ingredient is the use of a certain convenient gauge, space-cone. We have shown that each 
tree level gluon Feynman diagram “factorizes”, i.e. it can be written as a sum of product of 
lower on-shell amplitudes that arise from successive cuts of all internal lines which connect two 
reference gluons. To be able to still satisfy the momentum conservation laws, the momenta 
of the reference gluons must be shifted in a particular way. Moreover, the amplitudes arising 
from cuts must be multiplied by the propagator of the line which was cut. The proof of 
factorization is based, among other things, on an algebraic identity involving the momentum 
space propagators, which we recognized as the complexified Fourier transform of the largest 
time equation. 

Let us recall that the largest time equations are exact identities in quantum field theories, 
independent of loop levels. They yielded results for spectral representation of two-point func¬ 
tions, side-wise dispersion relations for vertex functions, etc. We may infer from these examples 
that exploiting them in the area of complexified unitarity to complement space-cone gauge free¬ 
dom should likely offer new opportunities at the loop level for organizing QCD computations, 
among other things. In particular, it should facilitate a field theoretical perspective of the 
one-loop recursion relations [11]. 

We have also addressed generalizations of the tree level BCFW recursion relations involving 
massive charged scalars and fermions, proving them in a purely field theoretical setup. 


Note added 

While this manuscript was being written, we became aware of [12] which has partial overlap 
with our results. 
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